Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 20 March 2013 (MN W^X style file v2.2) 



Flux calibration of the AAO/UKST SuperCOSMOS Ha Survey 

fr) David J. Frew 1 ' 2 *, Ivan S. Bojicic 1 ' 2 ' 3 , Quentin A. Parker 1 ' 2,3 , Mark J. Pierce 1 ' 4 , 
M.L.P. Gunawardhana 5 and WA. Reid 1,2 

1 Department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109 Australia 

2 Research Centre for Astronomy, Astrophysics and Astrophotonics, Macquarie University, Sydney, NSW 2109 Australia 
i 3 Australian Astronomical Observatory, PO Box 915, North Ryde, NSW 1670, Australia 

■ 4 Department of Physics, University of Bristol, Bristol, BS8 1TL, UK 

■ 5 Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006, Australia 



O 



Accepted . Received ; in original form 



ABSTRACT 

The AAO/UKST SuperCOSMOS Ha Survey (SHS) of the southern Galactic plane was, when 
completed in 2003, a powerful new addition to wide-field surveys. It has a combination of 
areal coverage, spatial resolution and flux sensitivity in a narrow imaging band which still 
marks it out today as an excellent resource for the astronomical community. The 233 separate 
fields are available online in digital form, with each field covering 25 square degrees. The 
SHS has been the motivation for equivalent surveys in the north, and new digital Ha surveys 
now beginning in the south such as VPHAS+. The SHS has been the foundation of many im- 
portant follow-up discovery projects in the southern sky with the Macquarie/A AO/Strasbourg 
Ha (MASH) planetary nebula project being a particularly successful example. However, the 
full astrophysical potential of the SHS has been hampered by lack of a clear route to accept- 
able flux calibration from the base photographic data. We have determined the flux calibration 
factors for 170 separate SHS fields, and present a direct pathway to the measurement of in- 
tegrated Ha fluxes and surface brightnesses for resolved nebulae detected in the SHS. We 
also include a catalogue of integrated Ha fluxes for 88 known PNe, and several Hii regions 
and other objects measured from the SHS, and use these data to show that integrated fluxes, 
accurate to ±0.1 - 0.15 dex, can be generally obtained across these survey fields. For the re- 
maining 63 fields, a mean calibration factor of 12.0 counts/px/rayleigh can be used, which 
will allow the determination of reasonable integrated fluxes accurate to better than ±0.2 dex. 
We outline the procedures involved and the problems and limitations that need to be appre- 
ciated in achieving such flux measurements. This paper forms a handy reference source that 
will significantly increase the scientific utility of the SHS. 

Key words: astronomical databases: surveys; techniques: photometric; ISM: general, plane- 
tary nebulae; HII regions; supernova remnants 



1 INTRODUCTION 

The AAO/UKST SuperCOSMOS Ha Survey (SHS; Parker et al. 
2005, hereafter Paper I) provides high-resolution (^1"), narrow- 
band, digital imaging data of the Southern Galactic Plane (SGP). 
The excellent combination of survey characteristics has made the 
SHS a powerful tool for the detection of extended, low surface 
brightness emission nebulosity over a range of angular scales from 
arcseconds to degrees (refer to later references). The SHS is also 
sensitive to faint point source Ha emitters such as WR stars, Be 
and B [e] stars, T Tauri stars, symbiotic stars, cataclysmic variables, 
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and other exotica (Drew et al. 2003; Hopewell et al. 2005; Pretorius 
& Knigge 2008). 

The survey utilised the 1.2-m UK Schmidt Telescope (UKST) 
at Siding Spring Observatory between 1998 and 2003. Exposures 
were taken through a high specification, monolithic Ha interfer- 
ence filter with a clear aperture of approximately 305 mm, a cen- 
tral wavelength of 6590A, a peak transmission of ^90% and a 
full width at half maximum (FWHM) of 70A (Parker & Bland- 
Hawthorn 1998). This was sufficient to cover the Ha and both [N n] 
AA6548,84 A emission lines across the full range of radial velocity 
observed for Galactic emission (±300 kms -1 ). 

The SHS footprint covers about 4000 square degrees of the 
SGP between I = 210°- 30°, and to a latitude, |b| «10 - 13°. It 
comprises 233 distinct but overlapping 5° survey exposures on non- 
standard field centres 4° apart (cf. the normal 5-degree spacing), 
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together with 40 fields covering a separate contiguous region of 
650 square degrees in and around the Magellanic Clouds (Morgan, 
Parker & Phillipps 1999). Each survey field consisted of a 3 -hour 
Ha Tech-Pan exposure and a matching, usually contemporaneous, 
15-min broad-band (5900-6900 A) 'short-red' (SR) exposure taken 
with an OG590 filtei[jfor continuum subtraction purposes. To date 
the SHS is the only high-resolution Ha survey covering the full 
SGP. Earlier work by J. Meaburn and co-workers used a square 
multi-element interference filter on the UKST, with a FWHM of 
105 A when used in combination with the 098-04 emulsion (El- 
liott & Meaburn 1976), but the plates only covered limited areas of 
the SGP (e.g. Elliott, Goudis & Meaburn 1976; Davies et al. 1978; 
Meaburn & White 1982). However, the coverage of the Magellanic 
Clouds was relatively complete (Davies, Elliott & Meaburn 1976; 
Goudis & Meaburn 1978; Meaburn 1980). However these survey 
plates are not generally available to the community. 

The scanned SHS data at 10/xm (0.67") resolution is avail- 
able online in digital form through the SHS website0, hosted by the 
Royal Observatory, Edinburgh. Note that the SHS coverage is over 
a factor of two larger than the equivalent northern INT Photometric 
Ha: Survey (IPHAS; Drew et al. 2005) which only extends to |b| 
~ 5°. There is a small area of overlap between these two surveys. 
The new VST/OmegaCam Photometric Ha Survey of the South- 
ern Galactic Plane (VPHAS+; Drew et al., in preparation jf| is the 
southern counterpart of IPHAS, and also only extends to |b| ~ 5°. 
VPHAS+ has been underway on the ESO VST telescope in Chile 
since the start of 2012, but will only supersede the utility of the 
SHS in the Galactic bulge and inner plane (|b| < 5°) regions. 

The SHS survey reaches a magnitude depth of approximately 
Rs — 20.5 for continuum point sources, and has a faint sensitivity 
limit to diffuse emission of 2-5 rayleighfl (see Paper I for fur- 
ther details). This can be compared to the approximate limit of the 
best Palomar Observatory Sky Survey (POSS) red plates, which 
used the 103a-E emulsion, of ~20 rayleighs (see Peimbert, Rayo 
& Torres-Peimbert 1975), a factor of 10 brighter than the SHS limit. 
Both the POSS-II (Reid et al. 1991) and UKST red plates used the 
IIIa-F emulsion; these plates have a sensitivity limit about twice as 
deep as the POSS, or ^10 - 12 R. In comparison, sky-limited nar- 
rowband CCD sensitivity limits are at the ^1 R level (e.g. Xilouris 
et al. 1996), comparable to the Southern Ha Sky Survey Atlas 
(SHASSA; Gaustad et al. 2001) and the Virginia-Tech Spectral-line 
Survey (VTSS; Dennison, Simonetti & Topasna 1998). For sum- 
maries of past, present and upcoming Ha surveys, see Paper I and 
Frew, Bojicic & Parker (2013, hereafter FBP13). 

An in-depth discussion of the technical details of the SHS was 
presented in Paper I. The scientific exploitation of the SHS is still 
underway, as it has proven to be an excellent discovery medium 
for both extended nebulae and compact Ha emitters. One area of 
particular success has been in the detection of unprecedented num- 
bers of new Galactic planetary nebulae (PNe), published in the two 
Macquarie/AAO/Strasbourg Ha (MASH) catalogues (Parker et al. 
2006; Miszalski et al. 2008). While large numbers of these PNe 
are concentrated in the Galactic bulge (Acker, Peyaud & Parker 



1 The effective wavelength of the short-red band is blueward of the widely 
used R F system (III-aF emulsion + RG630 filter; Couch & Newell 1980), 
hence the name. The colour equations of this combination have been out- 
lined by Morgan & Parker 2005. 

2 see http://www-wfau.roe.ac.uk/sss/halpha/ 

3 see http://www.vphasplus.org/ 

4 1 rayleigh = 10 6 /47r photons cm -2 s _1 sr _1 = 2.41 x 10 -7 erg cm -2 
s _1 sr _1 at Ha. 



2006), many are found throughout the SHS footprint along the SGP. 
These discoveries and their subsequent spectroscopic confirmation 
are having a major impact on PN science (see Parker et al. 2012a, 
for a review). 

The new MASH PNe, the recent IPHAS discoveries (e.g. 
Mampaso et al. 2006; Viironen et al. 2009, 2011; Sabin et al. 2010; 
Corradi et al. 2011), and other significant contributions (e.g. Ja- 
coby et al. 2010; Kronberger et al. 2012) have pushed the currently 
known Galactic population up to more than 3200 PNe (FBP13), 
more than double what it was a decade ago. The SHS also proved 
useful for finding previously uncatalogued emission features such 
as knots, jets and outer halos around known PNe (e.g. Miszalski 
et al. 2009; Frew, Bojicic & Parker 2012; Frew et al. 2013b) and 
symbiotic stars (e.g. Miszalski et al. 2012). 

The survey has likewise been a boon for uncovering more ex- 
tended nebulae such as the optical counterparts of known Galactic 
supernova remnants (e.g. Stupar, Parker & Filipovic 2007a; Stu- 
par et al. 2007b; Stupar & Parker 2009, 2011, 2012; Robbins et 
al. 2012) seen previously only in the radio domain. Painstaking 
searches of the SHS by Stupar, Parker & Filipovic (2008) also 
revealed many faint optically detected SNRs and SNR candidates 
not previously recognised in the radio (see also Walker, Zealey & 
Parker 2001; Alikakos et al. 2012). Similarly, the potential of the 
survey for deep optical detections of Hllregions and SNRs in the 
LMC has been realised by Muller & Parker (2007), Bojicic et al. 
(2007) and Reid & Parker (in preparation). 

Similarly, the SHS has also been used to study star-forming 
regions, both compact and diffuse (e.g. Cohen et al. 2002; Pierce 
2005; Cohen et al. 2007, 2011; Zavagno et al. 2007; Urquhart 
et al. 2007; Cappa et al. 2008; Deharveng et al. 2009; Pinheiro 
et al. 2010; Russeil et al. 2005, 2010; Bik et al. 2010; Murphy 
et al. 2010), the diffuse interstellar medium (ISM; Tackenberg et 
al. 2013), Herbig-Haro objects (Girart & Viti 2007), giant Herbig- 
Haro outflows (e.g. Mader et al. 1999), and a range of ejecta, shell 
and bowshock nebulae around massive stars (e.g. Cohen, Parker & 
Green 2005; Gvaramadze & Bomans 2008; Gvaramadze, Kniazev 
& Fabrika 2010; Stock & Barlow 2010; Stock, Barlow & Wesson 
201 1 ; Duronea, Arnal & Testori 2012; Duronea, Arnal & Bronfman 
2013). 

Finally, the SHS has been utilised for point- source studies re- 
vealing a range of emission line-stars including WR stars (Drew et 
al. 2003; Morgan, Parker & Cohen 2003; Hopewell et al. 2005), Be 
and T Tauri stars (Pierce 2005), and cataclysmic variables (e.g. Pre- 
torius & Knigge 2008). Importantly, using the related SHS imaging 
data for the LMC, Reid & Parker (2006a, b) used a deep stack of 12 
Ha and six matching SR exposures spanning the central 25 deg 2 of 
the LMC to uncover over 2,000 compact emission sources, includ- 
ing over 450 new PNe and a large number of emission-line and vari- 
able stars. These included WR stars, carbon stars, symbiotic stars, 
and other types of variable stars. The confirmed hot emission-line 
stars, 89% of which are type Be, were published together with their 
underlying spectral classification and luminosity type, radial and 
rotational velocities, Ha fluxes, equivalent widths and photometry 
by Reid & Parker (2012). 

1.1 Motivation for the Project 

Despite these studies, the full potential of the SHS has been ham- 
pered by the lack of a clear procedure to establish a reliable, inte- 
grated Ha flux or a surface brightness measurement for the range 
of distinct emission structures visible on the survey. Thus we need 
a reliable method to transform the pixel intensity values from Su- 
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perCOSMOS microdensitometer scans (e.g. Hambly et al. 1998) of 
the Ha images into intensity units such as rayleighs. The procedure 
should also be consistent from field to field. Pierce (2005) and Pa- 
per I showed that the survey data have been well exposed to capture 
Galactic emission on the linear part of the photographic character- 
istic curve (e.g. Smith & Hoag 1979; Parker & Malin 1999). How- 
ever, no detailed process for achieving reliable flux calibration for 
discrete sources on the SHS has yet been provided in the literature. 

To address this deficiency we extend this earlier work so re- 
searchers can derive meaningful results from flux-calibrated SHS 
data as required, following the precepts outlined in this paper. The 
motivation for this study is to explore the viability of using the 
SHS imaging data to establish a reliable mechanism for Ha flux 
measurements. This process is based on the continuum subtracted 
images of non-saturated, resolved emission sources not detected 
on the SHASSA (Gaustad et al. 2001), primarily due to the mod- 
est (48") resolution of the latter survey. This study also attempts 
to address the various limitations of the base photographic data. 
This includes the restricted dynamic range (when scanned by the 
SuperCOSMOS microdensitometer) compared to the inherent pho- 
tographic density range of the original films (e.g. MacGillivray & 
Stobie 1984; Hambly et al. 2001). In a few cases the fields are af- 
fected by high chemical fog (e.g. Tritton 1983) which further re- 
stricts the dynamic range of these survey exposures. 

The structure of this paper is as follows: §2 describes the lim- 
itations of the photographic survey data and §3 describes the flux 
calibration process for the SHS. In §4 we focus on the derivation 
of SHS Ha fluxes from the imaging data, and in §5 we undertake a 
comparison of the SHS fluxes for a range of discrete nebulae with 
independent measurements. Our conclusions and recommendations 
are presented in §6. 



2 FLUX CALIBRATION OF THE SHS SURVEY IMAGES 

The calibration of the SHS survey images is achieved by direct 
comparison of the continuum- subtracted SHS data with the match- 
ing continuum- subtracted SHASSA data (Gaustad et al. 2001) 
which provides narrowband (Ha + [Nil]) and continuum images 
of the entire southern sky. The preliminary flux calibration process 
was described by Pierce (2005) and Paper I, and the reader is re- 
ferred to the latter paper for technical details. Although SHASSA 
has a low spatial resolution with 48" pixels, it has the consider- 
able advantage of being continuum subtracted and accurately flux 
calibrated in rayleighs, using a range of calibrating PNe from Do- 
pita & Hua (1997). The SHASSA calibration has been verified by 
Finkbeiner (2003) who used the independent data from the Wis- 
consin H-Alpha Mapper (WHAM; Haffner et al. 2003) Northern 
Sky Survey at 1° resolution. Pierce et al. (2004), Pierce (2005), 
Frew (2008), and FBP13 also provided independent analyses of the 
SHASSA data. Indeed this work has recently culminated in the use 
of the SHASSA data to provide accurate, integrated Ha fluxes for 
over 1000 Galactic PNe (FBP13). 

In this paper we have completely redone the SHS calibration 
process afresh, by selecting data from the unvignetted ~25 square 
degree central regions (see Figure of all 233 SHS fields. This 
was undertaken using the onlin^fl 16 x blocked-down data for each 
field. These images have a resolution of 10.7" per pixel, and were 
not available in that format to Pierce (2005). These data provide a 

5 see http://www-wfau.roe.ac.uk/sss/halpha/hafields_online.html 



more convenient way of working across large areas for comparison 
with the calibrated SHASSA data. It is important to remember that 
one 13° diameter SHASSA field completely covers a single 5° di- 
ameter SHS field. Of course the SHASSA survey itself suffers from 
artefacts resulting from the effects of bright stars (see Gaustad et al. 
2001, and Figured). 

2.1 Continuum Subtraction 

Each matching SR image provide an effective measure of the con- 
tinuum component for each SHS field, despite including the Ha 
emission line in the bandpass. The Ha and SR images are gener- 
ally exposed to attain the same depth for continuum point sources 
(Paper I). However, the nature of photography and the vagaries of 
the observing conditions (e.g. if the exposure pairs were not con- 
temporaneous and taken in different moon phases or if the see- 
ing changed) mean that the SR and Ha point source depth im- 
age quality is similar, but not exactly so. Thus, it is necessary to 
precisely determine a continuum- subtraction scaling factor (SF) on 
a field-by-field basis for the continuum subtraction process to be 
effective. For high-dynamic range CCD exposures, the standard 
method for determining the appropriate scaling factor to subtract 
continuum from narrow-band is to compare aperture photometry 
for stars on both images whose exposures are normally interleaved 
on short time- scales. This does not work so well with the Ha plus 
SR film exposure pairs (Pierce 2005), often leading to under- or 
over- subtraction of the continuum. This produces positive or nega- 
tive residuals around point sources, so care is needed (see £14.21 ). 

Firstly, the 16 x blocked-down (10.7'Vpixel) Ha and SR sur- 
vey images for each SHS field were obtained from the SHS website. 
Then, a continuum- subtracted Ha(+ [N II]) image was created from 
the matching Ha and SR images in IRAF, using a scaling factor 
(SF), as defined in Pierce (2005) and Paper I. A range of values 
for the scaling factor can be applied to the SHS data until the zero- 
point of the continuum- subtracted SHS survey images best match 
the zero-point of the equivalent SHASSA data, indicating the ap- 
propriate value to use. The accuracy of the scaling factor is critical 
in order to avoid incorrect subtraction of the continuum compo- 
nent. Once a suitable SF has been determined, the continuum sub- 
traction process works well, removing most of the stellar images 
and the diffuse continuum. Only stars that sit in regions with the 
strongest Ha emission have been over- subtracted, and hence ap- 
pear as white spots in the image, as do any high-surface brightness 
knots in fainter extended nebulae. For a few fields, we adopt SF = 
1, if the comparison is inadequate to determine a formal value. 

For example, Figure Q] shows four images of a 80' square re- 
gion from field HAL0970, which shows Ha emission at a range 
of intensities. The left image is the 16 x blocked-down SHS Ha 
image, the second image is the continuum- subtracted Ha im- 
age, the third image shows this image re-binned to a resolution 
to match SHASSA, while the right image shows the equivalent 
SHASSA data. In this example, this continuum subtraction pro- 
cess has worked well, removing most of the stellar images and the 
diffuse continuum. Only the stars that sit in the strongest emission 
have been over- subtracted and appear as white spots in the image. 
Note the residuals around stars seen in the SHASSA image, and the 
planetary nebula M 1-46, marked with a circle. 

2.2 Calibration Factors 

Each SHS field has a corresponding calibration factor (CF) de- 
termined by Pierce (2005), but these are not generally available 
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Figure 1. Comparison of SHS and SHASSA data. From left to right: (I) Ha image of a 80' x 80' region near the centre of SHS field HAL0970; (2) matching 
continuum-subtracted Ha: image; (3) re-binned continuum- subtracted SHS Ha image; (4) equivalent flux-calibrated SHASSA Ha image. Note the stellar 
residuals in the SHASSA image, and the planetary nebula M 1-46, highlighted with a circle in each panel. A colour version of this figure is available in the 
online journal. 



to the community. Consequently, we have first re-calculated and 
improved the reliability of the CF for each survey field using our 
somewhat different approach based on the 16 x binned SHS data. 
In addition, in fields lacking obvious large scale emission, we have 
attempted to identify isolated nebulae (such as PNe) with which to 
provide a new calibration. 

For each SHS field, we selected the corresponding SHASSA 
field by minimising the distance between the respective field cen- 
tres. Secondly, the 16 x blocked-down continuum- subtracted im- 
ages were re-binned and re-aligned to match the correspond- 
ing 48 // resolution SHASSA image using the MONTAGE toolkit 
(http://montage.ipac.caltech.edu/). The SHASSA image was then 
cropped to the re-projected and re-binned SHS image. Now, 
the pixel intensity counts of the re-binned SHS images can be 
compared directly with the equivalent pixel intensity values (in 
rayleighs) from SHASSA to yield the CF for that field. The 
MONTAGE mPro ject module also returns the average value of 
binned pixels, which means that the final CF, derived from the 
pixel-to-pixel comparison with SHASSA, applies to the original 
full resolution SHS image. Large (typically 30' - 60') representa- 
tive areas from the normal unvignetted area of each SHS field were 
selected and matched to give a robust calibration for each field. 
We emphasise that these regions were generally available for cal- 
ibration, due to the large field-to-field overlap of the survey. This 
known vignetting from the Ha filter was the original motivation 
for the SHS field centres being placed on a 4° grid rather than the 
standard 5° grid, so this problem can be avoided (see Paper I). 

Figure |2] shows the blocked-down SHS Ha image of field 
HAL0970 with contours of filter response at 70%, 75%, 80%, 85%, 
90%, 92%, 95%, 97% and 98% of peak overlaid. Selected subre- 
gions are marked by rectangles of different colour. The magenta 
rectangles represent three subregions (numbered 1-3) inside the 
unvignetted part of the field, and the blue rectangle (region 4) in- 
cludes the bright Hll region M 17 (NGC 6618), also located within 
the unvignetted part of the field. For a comparison, the green, red 
and orange rectangles (regions 5-7) are areas of the field affected 
by filter vignetting. 

A graphical plot of continuum- subtracted SHS pixel intensi- 
ties against equivalent SHASSA pixel intensities should return a 
linear relation with a common zero point if the reduction and in- 
tensity calibration have been properly carried out and assuming 
the photographic intensity values are unsaturated and fall on the 
straight-line portion of the characteristic curve. The gradient of this 
plot is the CF for the image (see below). If the CF determined from 
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Figure 2. The full 5° x5° blocked-down SHS Ha image of field HAL0970 
with overlaid with contours of the filter response at 70%, 75%, 80%, 85%, 
90%, 92%, 95%, 97% and 98% of peak. Selected subregions are marked by 
rectangles as described in the text. Those numbered 1-4 are taken from the 
best area of filter response and were chosen to select a range of emission 
line intensities. Regions 5-7 probe the edges of the field where the filter 
transmittance drops below 85% and are not normally accessed by the online 
tool. A colour version of this figure is available in the online journal. 



the comparison plot is well behaved, it is applicable to full resolu- 
tion SHS data for the entire field. Unfortunately, some survey fields 
exhibited very little diffuse emission making calibration for these 
fields less straightforward. In such case a careful search for small 
areas of discrete, low intensity Ha emission in the unvignetted zone 
of the field was made, including any known large planetary nebulae. 
This enabled pixel by pixel comparisons between the two surveys 
to deliver appropriate CFs. 

The process is illustrated for field HAL0970 (Figure 0. We 
show the calibration plot of the selected rectangular subregions as 
Figure [3] The points from each subregion are plotted as a different 
colour. The magenta points are from three subregions (1-3) inside 
the unvignetted part of the field, and the blue points represent the 
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Figure 3. Calibration plot of the numbered and colour-coded subregions 
on field HAL0970. The magenta and blue points are from four subregions 
(1-4) inside the normal unvignetted part of the field, with the blue points 
representing the bright Hn region M 17, extending to higher intensities. 
These align with the magenta points on the one trend, but show the effect 
of bright photographic saturation, seen above 500 rayleighs. The green, red 
and orange points show data from areas of the field (regions 5-7) affected by 
vignetting, and these are displaced from the magenta/blue points. A linear 
fit to the unvignetted data (between 50 and 500 R) is also plotted (red dashed 
line), with a gradient (CF) of 14.0 SHS pixel counts per rayleigh. A colour 
version of this figure is available in the online journal. 



Figure 4. Calibration plot of field HALO 175 which is a field of intense Hen 
emission centred on the Eta Carina Nebula. Owing to this the Ha exposure 
was reduced to two hours. The plot illustrates beautifully the effect of flux 
oversubtraction at high surface brightness. A linear fit to the data between 
50 and 700 R is also plotted, with a gradient (CF) of 1 1.2 SHS pixel counts 
per rayleigh. A colour version of this figure is available in the online journal. 



bright Hll region M 17 (region 4), also located within the unvi- 
gnetted part of the field. These align with the magenta points on the 
one trend, but show the effect of bright photographic saturation, 
seen above 500 rayleighs. A least-squares linear fit to the magenta 
and blue data points below 500 R is also plotted, with a gradient of 
14.0 SHS pixel counts per rayleigh. This demonstrates that if pro- 
vided areas are chosen from within the unvignetted field, consistent 
CFs are obtained giving confidence in its applicability across the 
field. However, if pixel-pixel comparisons are made from the sub- 
regions in Figure |3] which sample data from the vignetted regions 
of the field (green, red and orange points; regions 5-7), the agree- 
ment is poor and varied. The spread in SHS counts is artificial and 
results from applying the flat field correction to the underexposed 
areas of film. 

We also illustrate the calibration plot for HALO 175 which is 
a field of intense Ha emission centred on the Eta Carina nebula 
(Figure |4}. Because of this, the Ha exposure was reduced to two 
hours, so photographic saturation sets in at a higher than usual level 
of 700 - 800 R. The plot illustrates beautifully the effect of flux 
oversubtraction at high surface brightness, as the SR 'continuum' 
image also saturates in the brightest section of the Carina nebula. A 
linear fit to the data between 50 and 700 R is plotted as a red dashed 
line, with a gradient (CF) of 1 1.2 SHS pixel counts per rayleigh. 

This calibration process was repeated for all 233 SHS fields. 
As a result, more fields now have reliable CFs than were deter- 
mined by Pierce (2005) and the data are presented in Table [T] The 
comparison plots for six well-behaved SHS fields are shown in Fig- 
ure [5] Here an essentially linear pixel-intensity relation can be seen 
between the SHS and SHASSA data and the linear fits are overlaid. 
A further six fields with some caveats are illustrated in Figure [6] 

In Figure [6] we show the comparison plots for six second- 



ranked SHS fields, which have linear fits that are not as precise as 
the first-ranked fields. In the case of HAL1018, subtle changes in 
the [N Il]/Ha ratio across an ionization front in the large Hll region 
IC 2177 produce a variable ratio of SHS pixel counts / R across the 
nebula. Furthermore, in other fields, discrete but widely separated 
Hll regions may occur in areas of the filter response which vary at 
the 2-3% level. HALO 137 is such a case, where the slight incon- 
sistencies around the linear fit seen are due to the slightly varying 
filter response in different areas of the field where the discrete Hll 
regions occur. HAL1 150 is an example of a field with a limited in- 
tensity range of emission, while the other fields show varying levels 
of photographic fog. These various effects can conspire to produce 
a formal fit with apparent curvature, or the appearance of two or 
more closely spaced trends on the calibration plot. These effects 
are usually small however, and should not compromise the utility 
of the fit in these SHS fields. The CFs determined from the linear 
fits can be applied with some confidence to the full (0.67") resolu- 
tion pixel data for the relevant survey field. 

The SHS coverage is illustrated in Figure |7] where the quality 
rank for each field is colour-coded. Note that the third-ranked fields 
are essentially devoid of diffuse Ha emission. Figure [8] shows the 
CF distribution for 170 fields (of 233) that could be fit with a linear 
function. The distribution of the first ranked-fields has a mean value 
of 11.9 ± 2.9 counts/px/R (la dispersion), while the second-ranked 
fields have a mean 12.2 ±3.2 counts/px/R. The weighted mean is 
12.0 ± 3.0 counts/px/R, which can be used as a default value for 
those fields without a formal CF determination. This approach then 
provides a direct mechanism to determine Ha fluxes for individual 
emission structures not resolved by SHASSA, including, for exam- 
ple, many new compact MASH PNe discovered from the SHS data 
(e.g. Parker et al. 2006; Miszalski et al. 2008). 
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Table 1. Data summary for all 233 SHS survey fields. Note there are 63 fields with no formal CF determined (usually due to poor dynamic range or high 
chemical fog, or both). A default value of 12.0 ± 3.0 pixel counts / R can be used for these fields, as explained in the text. The table is published in its entirety 
as an online supplement, and a portion is shown here for guidance regarding its form and content. 



Field 


SF 


CF 

(count/px/R) 


Quality 


Comment 


HAL0067 




(12.0) 


3 


small dynamic range 


HAL0068 




(12.0) 


3 


small dynamic range 


HAL0069 




(12.0) 


3 


small dynamic range 


HAL0070 




(12.0) 


3 


photographic fog 


HAL0096 




13.3 


3 


small dynamic range 


HAL0097 


0.98 


11.4 


2 




HAL0098 


0.83 


13.9 


2 


photographic fog 


HAL0099 


0.42 


(12.0) 


3 


photographic fog 


HAL0100 


1.10 


12.7 


1 




HAL0101 


1.12 


11.7 


2 


small dynamic range 


HAL0102 




(12.0) 


3 


small dynamic range 


HAL0103 


0.44 


6.2 


2 


photographic fog 


HAL0131 


0.62 


(12.0) 


3 


photographic fog 


HALO 132 


0.91 


10.0 


2 




HAL0133 


1.29 


6.3 


2 




HALO 134 


1.27 


10.3 


1 


saturates at 500 R 


HALO 135 


0.92 


6.6 


1 




HAL0136 


0.87 


12.0 


2 




HALO 137 


1.24 


18.9 


2 


large scatter 


HAL0138 


0.91 


8.7 


2 


large scatter 
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Figure 5. SHASSA-SHS pixel-by-pixel comparison plots for six first-ranked SHS fields. Some SHS fields (e.g. HAL0350, HAL0482 and HAL0555) show 
evidence of saturation where the intensity exceeds ~500-600R. Below this level, most fields show a linear relation between the SHASSA and SHS. Note 
that fields HALO 134 and HAL 1060 show evidence of photographic fog, which elevates the sky level on these SHS fields. A colour version of this figure is 
available in the online journal. 



3 THE LIMITATIONS OF THE SHS SURVEY DATA 

As previously mentioned not all SHS fields provide reliable flux 
measurements Some fields are also naturally devoid of extensive 
Ha emission, especially those at higher Galactic latitudes. In these 
cases it can be problematic to find unvignetted areas with emis- 



sion features that demonstrate sufficient dynamic range to provide 
a large enough lever-arm for a reliable calibration. 



Unlike CCD data, photographic exposures have a limited lin- 
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HAL630 



500 600 




HAL1150 



SHASSAIR] SHASSA[R] SHASSA [R] 

Figure 6. SHASSA-SHS pixel-by-pixel comparison plots for six second-ranked SHS fields. The slight inconsistencies around the linear fit seen in HAL0137 
are due to the slightly varying filter response in different areas of the field where the discrete Hn regions occur. HAL 1150 is an example of a field with a 
limited intensity range of emission, while the other five fields show varying levels of photographic fog. 




22 20 18 16 14 12 10 8 6 

RA (J2000) 



Figure 7. The sky coverage of the SHS in equatorial J2000 coordinates (note that truncated RA and DEC coordinates). The green, amber, and white boxes 
represent first-, second-, and third-ranked fields, respectively, and the field numbers are given. A 'clickable' version (though without colour-coding) is available 
on the SHS website. A colour version of this figure is available in the online journal. 
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Calibration Factor 



survey fields where each SHS field was ranked according to suit- 
ability for obtaining accurate Ha fluxes. These ranks are given in 
Table Q] Representative samples of first- and second-ranked fields 
are shown in Figure [5] and Figure [6] 

In summary, of the 233 SHS fields, 50 (21%) are well con- 
strained by a linear fit below 500 R and so have the top ranking of 
1, while a further 120 fields (52%) also have a linear fit, but with 
some caveats, and have been given a quality rank of 2. The 63 re- 
maining SHS fields (27%) have little diffuse emission, making the 
determination of reliable CFs less straightforward, primarily due to 
a very restricted dynamical range of emission. A default value of 
12.0 pixel counts / R can be safely used for these fields. For ease of 
use, we give the CF derived from the linear fit for 170 SHS fields 
in Tabled] 



Figure 8. Histogram of derived calibration factors (CFs) for 170 SHS fields. 
The green and orange bars represent the numbers of first- and -second- 
ranked fields, respectively. A colour version of this figure is available in 
the online journal. 

ear response and saturate at a photographic densit)|f| of D ~ 4 (e.g. 
Hambly et al. 2001). However, the automated SuperCOSMOS mi- 
crodensitometer used to digitise the photographic exposures im- 
poses its own limitations, restricting the density range even further 
to D ^ 2 (Hambly et al. 1998). The data scanning and calibra- 
tion processes are affected not only by the relatively low dynamic 
range of SuperCOSMOS, but also because some Ha and SR pairs 
are not contemporaneous, being taken on different nights, and un- 
der different seeing conditions and lunar phase. Furthermore, these 
differences in exposure date may mean different hypersensitised 
emulsion batches were used which may have different levels of 
chemical fog and speed. Some fields have a higher than average 
level of photographic chemical fog due to the Tech Pan films be- 
ing push-processed. These factors can introduce variations in back- 
ground and sensitivity from exposure to exposure within the same 
'matching' Ha and SR pair. 

Further uncertainties in the emission levels across survey 
fields can arise from the presence of continuum emission, night-sky 
auroral lines, and geocoronal Ha emission as well as the effects of 
variable atmospheric extinction (Paper I). Ideally, all of these com- 
ponents would need to be disentangled to extract the true Galactic 
Ha emission. In practice, the geocoronal and auroral emission is 
considered as a low-level, temporally varying wash, which elevates 
the general background on each exposure to a slightly varying de- 
gree. This is considered to be less than 2 rayleighs across the SHS 
(Nossal et al. (2001), based on Ha data from the WHAM Northern 
Sky Survey (Haffner et al. 2003). These observations are being ex- 
tended to the southern sky as the WHAM-South Survey (Haffner et 
al. 2010). 

The original set of SHASSA and re-binned, continuum- 
subtracted SHS pixel-by-pixel comparison plots for selected SHS 
fields given by Pierce (2005) demonstrated the spread of dynamic 
range and linearity exhibited by the SHS data on a field by field ba- 
sis. Many fields have a well defined linear trend but some are less 
well-behaved due to some of the factors previously mentioned. As a 
consequence fields were allocated a quality flag (described below) 
to reflect this diversity. This is based on a fresh examination of all 

6 The photographic density, D, is the logarithm of the ratio of incident light 
to transmitted light through a region of a developed emulsion. 



4 RELIABILITY OF THE FLUX CALIBRATION 

In the previous section we have established revised calibration fac- 
tors for each SHS field. In order to demonstrate the capability 
of these fields to provide accurate integrated fluxes for discrete 
sources such as PNe, Hll regions, WR shells, and Herbig-Haro ob- 
jects, we identified a variety of suitable calibration objects across 
a range of SHS fields that have independent integrated Ha fluxes 
from the literature, or determined by us from SHASSA measure- 
ments. This process is described below. 

4.1 Selection of Individual Sources 

The object sample chosen for this flux calibration study includes a 
variety of astrophysical sources exhibiting a range of angular scales 
and emission intensities across the full range of SHS field expo- 
sure characteristics. These should be sufficient to determine if the 
calibration factors derived in the previous section are correct. The 
selected objects cover a Ha flux range of over ^7dex, and range 
in angular size from a few arcseconds to nearly 2° across. 

The sample comprises ^80 PNe, several Hll regions and 
Wolf-Rayet shells, and an isolated filament of the Vela supernova 
remnant. For each object the SHS Ha flux is measured follow- 
ing the recipe described in §4 and is either compared with the 
SHASSA-derived Ha flux from FBP13, or other independent Ha 
fluxes from the literature where available. For the ejecta nebulae 
Hf 39 (Hen 3-519) and Wray 15-751, the Ha fluxes, uncorrected 
for reddening, have been derived from the reddening-corrected 
fluxes and extinction values given in Stahl (1986, 1987), Hutse- 
mekers & van Drom (1991) and Hutsemekers (1994). 

4.2 Measurement Process 

Suitably sized Ha and matching SR areas of each SHS field that 
contain the selected calibrating objects were either downloaded as 
full 0.67 // /pixel FITS images from the SHS website or extracted 
from the 16 x blocked-down data depending on the object's angular 
size. The SR image was then subtracted from the Ha equivalent in 
I RAF using the scaling factor provided in Table Q] This is straight- 
forward thanks to the in-built World Coordinate System (WCS) for 
each field. These continuum- subtracted Ha images were then used 
for flux measurements. 

A photometric aperture was defined to carefully enclose 
the entire object in I RAF, and sky background subtraction was 
achieved by using a suitable sky annulus positioned to avoid stellar 
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residuals and other artefacts. Unlike SHASSA, the superior reso- 
lution of the full-resolution SHS data often enables the fine struc- 
tures and extended halos of most target sources to be seen. The 
background-corrected photometry counts from the aperture were 
then used to measure the source Ha flux using algorithms similar to 
those described in FBP13. In Figure [9] we present two examples to 
illustrate the fitting procedure we employ to derive sky- subtracted 
integrated SHS counts for our calibrating sources. 

A problem can be the presence of stellar residuals within the 
image area of extended objects. In the Galactic plane, the stellar 
number density is high so residuals may also be present in the im- 
mediate region surrounding the chosen calibrator (Figure [TO]). Stel- 
lar residuals can sometimes appear around bright, saturated stars 
or compact emission sources in the continuum subtracted products 
(Pierce 2005; Paper I). Other challenges include residuals of the 
same angular size or larger than some compact PNe which partly 
overlap the PN image. Although these saturated residuals were re- 
moved via the same method as above, there can be some contami- 
nation of the PN emission. 

Although, stellar residuals were present in the images of most 
objects, not all intensity contributions from these were corrected. 
This was because the majority of the calibrating objects studied 
(mainly PNe) were of moderate to large angular extent and conse- 
quently their images were populated by many small stellar residu- 
als. Correcting for such large numbers of stellar residuals is difficult 
and time-consuming, and in many cases the effect of these residuals 
are accounted for in the sky subtraction process. 

Furthermore, the SR exposures also record Ha emission, 
which, when subtracted from the matching Ha exposures, can lead 
to an under estimation of the true Ha flux. In most cases this pro- 
cess appears to be self-compensating for the addition of Ha emis- 
sion by residuals, as our calibration results testify (see below). The 
integrated fluxes seem little affected for nebulae below a intensity 
limit of ^7000 counts per pixel. If the contribution to the total Ha 
emission of a source from stellar residuals seems significant, we 
have attempted to apply a correction so as to provide the most ac- 
curate possible continuum subtracted Ha fluxes. A correction was 
also applied due to the presence of a particularly bright central star, 
or if the calibrating object is located close to a bright, saturated 
star or stars being embedded in the object of interest. In the case 
of Hf 39, the ionizing star is prominent relative to the nebula (Fig- 
ure [TO]). We simply measured the flux of the star with a small aper- 
ture and subtracted this from the total flux. 



4. 2. 1 Deconvolution of [Nil] Emission 

Another factor that must be considered in determining Ha fluxes 
from narrowband imaging is the close proximity of the [Nil] 
A6548A and A6584A emission lines either side of Ha which are 
included in the SHS filter bandpass (Paper I). This is the case for 
most Galactic sources of modest radial velocity. The contribution 
from the [N II] lines to the total signal through the filter can vary 
significantly among emission sources. For example, for Type I PNe 
(Peimbert 1978), the [N Il]/Ha ratio is usually more than three and 
can sometimes exceed 10 (Ruiz 1983; Corradi et al. 1997; Perinotto 
& Corradi 1998; Kerber et al. 1998; Tajitsu et al. 1999; Frew, Parker 
& Russeil 2006). This must be taken into account when comparing 
SHS Ha fluxes with those determined independently in the litera- 
ture. This can be easily done when slit spectroscopy is available as 
this provides the observed [N ll]/Ha ratio. We have chosen calibra- 
tors that have spectra available, so the [N ll]/Ha ratios are known. 
These are taken from the literature sources listed in FB PI 3, or else- 



where. The observed ratios, usually obtained from long-slit spec- 
troscopy, are deemed representative of the whole object which is a 
safe assumption in most cases, though significant internal variation 
across some PNe and SNRs has been observed. 

The calculation of the SHS filter transmission properties is 
made more difficult by the blue shifting of the band pass with re- 
spect to incident angle in the fast f/2.48 converging beam of the 
UKST (see Parker & Bland-Hawthorn 1998). The effect was mod- 
elled by breaking down the contribution from the converging light 
cone into a series of concentric rings as described in Paper I. The 
shifting responses from each ring of the beam are weighted accord- 
ing to the areas of the ring and summed to generate a smeared filter 
response. This allows the calculation of the contributions from the 
[N II] lines to the total flux recorded by the survey. The reader is 
referred to Figure 3 of Paper I for an illustration of the modelled 
Ha filter transmission profile. Based on the smeared-out profile, 
the SHS filter transmits Ha at 80% and [Nil] A6548 and A6584A 
at 82% and 50% respectively. Given that the intrinsic ratio between 
the A6548A and A6584A [N II] line strengths is quantum mechani- 
cally fixed at ^3.0, the SHS filter transmits approximately 58% of 
[Nil], compared to 80% transmission of the Ha flux. This gives a 
[N Il]/Ha filter transmittance ratio of 0.725 which is the constant in 
equation [2 below. 

While the filter transmission falls towards the off-axis edge of 
each survey field (see Figure|2j, this is not considered significant in 
this study due to the purposely large overlap between adjacent SHS 
fields on a 4-degree grid (Paper I). This enables the data of interest 
to be preferentially selected from the non- vignetted central region 
of an available field. The filter transmission constant will vary at 
large off-axis angular distances but this is only relevant for the few 
fields at the eastern edges of the SHS region. 



4.3 SHS Ha Fluxes 

The equation used to derive the integrated SHS Ha+[N II] (or 'red') 
flux is essentially the same as the expression given in FBP13. The 
SHS red flux in cgs units is given by: 



F red = 5.66 x 10" 18 x Ps x 



/SUM\ 

I^fJ ergcm 



(1) 



where the first constant in the expression is the conversion 
factor from rayleighs to cgs units, Ps is the 'plate scale' of the 
digital SHS data (0.67"per pixel for the full resolution images, 
and 10.7 /x for the 16 x blocked-down images). The summed counts 
obtained from the aperture photometry procedure (SUM) need to 
be divided by the relevant calibration factor, CF, taken from ta- 
ble [TJ in order to determine the integrated red flux in units of 
ergcm _2 s _1 . The SHS field number can be obtained from the 
header of any downloaded FITS image. 

The integrated Ha flux is then obtained by correcting for the 
contribution from the contaminant [N II] lines (Frew 2008; FBP13) 
and is given by the following expression: 



F(Ha) S HS 



Ktr R[N II] + 1 



(2) 



where i?[Nii] is the observed spectroscopically determined 
[Nll]/Ha ratio for the emission object in question (mainly PNe) 
and the constant, K tT = 0.725, takes into account the throughput of 
the SHS filter. Following FBP13, the [N II] flux refers to the sum of 
the A6548 and A6584 lines. If only the brighter A6584 line is mea- 
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Figure 9. SHS continuum- subtracted images showing the photometric apertures (red circles) and background sky annuli (yellow circles) for three PNe: 
IC 1295 (left), SB 4 (middle) and the compact object JaSt 36 (right). The images are 400", 135", and 50" on a side, respectively, and have NE at top left. A 
colour version of this figure is available in the online journal. 




Figure 10. SHS continuum- subtracted images illustrating the effect of a bright stellar residual from the massive WN11 star, Hen 3-519 (in its surrounding 
nebula, Hf 39, left), several residuals, including from the binary central star of the toroidal PN, PHR J1818-1526 (middle), and a very bright, saturated PN with 
flux over-subtraction in its interior (NGC 2440, right). The catalogued diameter is given by the blue dashed circle in each figure, and the photometric apertures 
and background sky annuli are selected as before. The images are 170", 320", and 230" on a side, respectively, and have NE at top left. A colour version of 
this figure is available in the online journal. 



surable in a spectrum, F[Nll] is estimated as 1.33 x F(A6584) as 
the 6584 to 6548 A[N II] line ratio is 3:1. 

The uncertainty budget for the SHS Ha flux is determined fol- 
lowing the procedure presented in FBP13. Additionally, the uncer- 
tainty in the CF for the particular SHS field, as measured from the 
fit, was added in quadrature to the uncertainty in aperture intensity 
counts. 

We emphasise that the [N II] transmission differs between the 
SHS and the SHASSA survey filters (70A and 32A FWHM re- 
spectively), so the SHS filter passes considerably more [N II] emis- 
sion. The transmissible [N Il]/Ha ratios of the two survey filters are 
0.725 and 0.375 respectively. In the original calibration of the SHS 
from SHASSA data (Paper I), both Hll regions and the diffuse ion- 
ized gas were assumed to have negligible [N II] emission. However, 
the median [N ll]/Ha ratio for a large number of Galactic Hll re- 
gions is actually ^0.4 (derived from the data graphically presented 
by Frew & Parker 2010), which is in reasonable agreement with 
the average A6584/Ha ratio of ~0.3 seen in a sample of faint Hll 
regions (Madsen et al. 2006a). These authors also noted that the 
A6584/Ha ratio in the warm interstellar medium (WIM) generally 
increased below an Ha intensity of 10 rayleighs, approaching unity 
at the faintest measured intensities. Owing to these points, there is 
an offset of ^10-15 per cent in the flux obtained from using the CF 



values presented in Table Q] We therefore recommend that the for- 
mally determined SHS Ha fluxes given by equation |2]be increased 
by this factor; i.e. the logarithmic Ha fluxes should be brightened 
by 0.05 dex. 

The total uncertainty on the SHS Ha flux is given by the fol- 
lowing expression, which includes a SHASSA uncertainty term as 
the SHS is bootstrapped to this survey: 

2 2 2 2 2 

^(Ha) — Cphot + &R + 0"cal + ^SHASSA (3) 

where a p hot is the proportional uncertainty in the photometry 
measurement, or is the uncertainty in the deconvolved Ha flux 
due to the uncertainty in „] , cr ca i is the calibration uncertainty 
for the respective SHS field, and ctshassa is the known SHASSA 
calibration uncertainty of 8% (Gaustad et al. 2001). 



4. 3. 1 Emission Measures 

The SHS pixel data can also be used in the study of the warm inter- 
stellar medium (WIM) and low-surface brightness Hll regions on 
spatial scales down to 3-5 showing fine detail that is not seen on 
SHASSA. For ionized hydrogen, the surface brightness, Sua (in 
rayleighs) and the emission measure (EM, in units of cm _6 pc) are 
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related by the following expression (Haffner et al. 2003; Madsen, 
Reynolds & Haffner 2006a; cf. Reynolds & Ogden 1982): 

EM = J n e dl = 2.75 1 9 Sh<* e 2 - 2E{B ~ v) cm~ 6 pc (4) 

where t is the electron temperature in units of 10 4 K. As be- 
fore, the pixel counts can be converted into rayleighs using the CF 
for the specific SHS field, and then the contribution from the con- 
taminant [N II] lines need to be deconvolved from the total mea- 
sured intensity to determine an Ha surface brightness or emission 
measure (see £14.2.11 ). The resulting data has wide application to the 
study of the structure and ionization of Hi I regions and the WIM. 



5 COMPARISON OF SHS Ha FLUXES AGAINST 
INDEPENDENT ESTIMATES 

Table |2] gives our measured SHS Ha fluxes for selected PNe and 
other calibrating emission-line sources. The fluxes are based on our 
newly determined individual field calibrations. We also list the in- 
dependent integrated Ha fluxes for these objects. The SHASSA 
Ha fluxes are primarily taken from FBP13, or measured here in 
the same manner, while Ha fluxes for the other nebulae are drawn 
from Stahl (1987), Hutsemekers (1994), Nota et al. (1995), Duer- 
beck & Benetti (1996), Weinberger, Kerber & Grobner (1997), Pol- 
lacco (1999), Beaulieu, Dopita & Freeman (1999), Boumis et al. 
(2003, 2006), Jacoby & Van de Steene (2004), Ruffle et al. (2004), 
Madsen et al. (2006b), Parker et al. (2011, 2012b), and Frew et al. 
(2013b). FBP13 showed that the flux uncertainties from Cappel- 
laro et al. (2001) are substantial, so we have not used these in our 
comparison. 

FigurefTTIpresents our SHS Ha fluxes plotted against indepen- 
dent measurements from the literature for all 101 objects in the cal- 
ibrating sample. A tight linear relation is obtained over a range of 
7 dex (10 million) in integrated Ha flux. The lower panel plots the 
offset (in the sense of SHS minus literature) against the SHS flux. 
Unsurprisingly the first-ranked fields permit the measurement of 
the most accurate fluxes, and the measurement accuracy is slightly 
less for second ranked fields. Our fluxes show a dispersion with 
the literature data of 0.12 dex from first-ranked fields, and 0.15 dex 
for second-ranked fields. Considering the uncertainties on the lit- 
erature fluxes, this translates to a nominal uncertainty of 0.10 dex 
and 0.14 dex, respectively. This shows that photographic data, if 
carefully considered is a useful source for emission-line intensities 
and fluxes (see also the flux data from Abell (1966), as reduced by 
FBP13). 

In Figure Q21 we show the offset between the SHS and liter- 
ature Ha flux plotted as a function of Ha surface brightness (in 
erg cm -2 s _1 sr _1 ), in turn derived from the literature flux and the 
angular size of the object. The sample is comprised of all the ob- 
jects from Table |3 plus another 45 saturated PNe, including very 
bright PNe like NGC 3918 and NGC 2440 (Figure[L0). The fluxes 
for the bright PNe are taken from FBP13 (and references therein) 
and the angular diameters are from Tylenda et al. (2003), Frew 
(2008), or the references given in the footnotes to Table The 
twin effects of photographic saturation and flux over- subtraction 
occur for the brighter nebulae with a surface brightness, logS'(Ha) 
^ —4.0 erg cm -2 s _1 sr _1 , where an inflection pont in the plot is 
clearly seen. 

In a follow-up paper (Bojicic et al., in preparation), we will 
provide new Ha fluxes for all measurable PNe found in the SHS 




-l.o L 



-16 -15 -14 -13 -12 -11 -10 -9 -8 -7 

log F(Ha) SHS 

Figure 11. SHS versus literature Ha fluxes for the whole sample of un- 
saturated calibrating sources. Fluxes derived from first- and second-ranked 
fields are plotted with filled and open circles, respectively. The agreement 
is excellent over a range in flux of more than 10 million. The lower panel 
plots the delta flux, in the sense of SHS minus literature fluxes. 



footprint, that cannot be measured from SHASSA due to their faint- 
ness and/or compact size. These fluxes will complement the recent 
flux catalogue of FB PI 3. In particular, the new SHS Ha fluxes for 
faint Galactic bulge PNe will be valuable, as they can be compared 
with new radio continuum fluxes (e.g. Bojicic et al. 201 la,b) to esti- 
mate extinction values, or compared with integrated fluxes in other 
wavebands; e.g. the [O III] fluxes of Kovacevic et al. (2011), in or- 
der to investigate the faint end of the bulge PN luminosity function 
in different emission lines (e.g. Ciardullo 2010). 



6 POINT SOURCE PHOTOMETRY 

For completeness, we give a brief account of the SHS as a provider 
of photometric data for point- source emitters. Each Ha and match- 
ing SR survey field has been automatically processed through 
the standard SuperCOSMOS photometric and image processing 
pipelines (e.g. Hambly et al. 1998, 2001). These provide so-called 
IAM (image-analysis mode) data for detected point-source and 
compact extended sources in each survey passband above a locally 
thresholded background. An image-deblending algorithm based on 
areal profile fitting at image levels equally spaced in log (intensity) 
is also applied (e.g. Beard, McGillivray & Thanisch 1990; Ham- 
bly, Irwin & MacGillivray 2001) which can be useful in crowded 
fields such as the Galactic plane. Consequently, large catalogues 
of detected sources, including accurate positions, photometry and 
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Table 2. A comparison of our SHS Ha fluxes for 88 PNe and 13 other nebulae with independent integrated Hen fluxes taken from the literature, as explained 
in the text. The table is published in its entirety as an online supplement, and a portion is shown here for guidance regarding its form and content. 

Name logF(Ha) Reference logF(Ha) SHS Field Remarks 

(other) (SHS) 



Abell 18 


-12.00 ±0.09 


FBP13 


-11.92 ±0.09 


HAL 11 97 


Abell 23 


-11.86 ±0.09 


FBP13 


-11.82 ±0.08 


HAL0526 


Abell 26 


-12.24 ±0.08 


FBP13 


-12.31 ±0.15 


HAL0601 


Abell 27 


-12.01 ±0.05 


FBP13 


-12.01 ±0.15 


HAL0602 


Abell 44 


-11.49 ±0.10 


FBP13 


-11.56 ±0.08 


HAL0970 


Abell 45 


-11.26 ±0.10 


FBP13 


-11.35 ±0.15 


HAL1060 


Abell 48 


-11.58 ±0.06 


FBP13b 


-11.55 ±0.10 


HAL1241 


Abell 58 


-12.28 ±0.06 


FBP13 


-12.32 ±0.11 


HAL 1333 


BMP J1613-5406 


-11.60 ±0.09 


FBP13 


-11.50 ±0.12 


HAL0289 


BMP Jl 808- 1406 


-11.68 ±0.09 


FBP13 


-11.78 ±0.14 


HAL0969 


CGMW 3-2111 


-12.38 ±0.10 


B06, FBP13 


-12.40 ±0.10 


HAL1058 


CVMP 1 


-11.50 ±0.10 


FBP13a 


-11.63 ±0.14 


HAL0232 


DPV 1 


-12.5 ±0.3 


DB96, P99 


-12.50 ±0.11 


HAL0968 


FP J071 1-2531 


-10.69 ±0.04 


FBP12 


-10.51 ±0.11 


HAL0756 


FP J1824-0319 


-10.40 ±0.10 


M06, FBP13 


-10.44 ±0.10 


HAL1240 


Fr2-7 


-10.82 ±0.10 


FBP13 


-10.86 ±0.12 


HAL0097 


G4.4+6.4 


-11.25 ±0.06 


FBP13 


-11.31 ±0.10 


HAL0878 


GL1PN J1530-5557 


-13.86 ±0.10 


PC12 a 


-14.17 ±0.12 


HAL0233 


GL1PN J1557-5430 


-14.79 ±0.15 


PC12 a 


-14.70 ±0.15 


HAL0234 


GL1PN J1642-4453 


-14.94 ±0.15 


PC12 a 


-14.88 ±0.15 


HAL0413 



Sakurai's object 



Reference codes: ASTR91 - Acker et al. (1991); B03, B06 - Boumis et al. (2003, 2006); BDF99 - Beaulieu et al. (1999); DB96 - Duerbeck & Benetti 
(1996); F13 - this work; FBP13 - Frew et al. (2013a); FBP13b - Frew et al. (2013b); H94 - Hutsemekers (1994); HDM98 - Hua, Dopita & Martinis (1998); 
JaSt04 - Jacoby & Van de Steene (2004); KH93 - Kistiakowsky & Helfand (1993); FPR06 - Frew et al. (2006); M06 - Madsen et al. (2006); P99 - Pollacco 
(1999); PC12 - Parker et al. (2012); PFM1 1 - Parker et al. (201 1); RZ04 - Ruffle et al. (2004); S87 - Stahl (1987); SK89 - Shaw & Kaler (1989); WKG97 - 
Weinberger et al. (1997); XP94 - Xilouris et al. (1994); XPPT - Xilouris et al. (1996). 
Notes: a Fluxes revised using updated calibration factors. 

other image parameters are provided for each field from the SHS 
websit^ll A significant advantage of the SHS data is the ability to 
detect point sources that have been photometrically calibrated to 
CCD standards (e.g. Boyle, Shanks & Croom 1995; Croom et al. 
1999). The narrowband Ha images are also calibrated to an 'R- 
equivalent' scale. 

The 3-hour Ha and 15-min broad-band SR exposure times 
are chosen so that both reach similar depths of Rs — 20.5 for point 
sources (Paper I). Where an object is detected in one band but not in 
the other, a default value of 99.999 is given in the online IAM cat- 
alogue data for the magnitude in the missing bandpass. Positional 
and magnitude-dependent errors are corrected for in the data that 
are available through the SHS website. Photometric consistency is 
achieved by using the field overlap regions to establish zero-points 
across the survey. These corrected magnitudes provide a means, 
for example, of selecting point- source Ha emitters by looking for 
significant brightening on the Ha source compare to its SR coun- 
terpart, especially when combined with the SSS /-band photome- 
try which effectively eliminates late-type stellar contaminants (Pa- 
per I). The IAM data was successfully used by Pierce (2005) to 
discover a range of emission-line stars, by Hopewell et al. (2005) 
who found several WO and WC stars, and in the discovery of nu- 
merous compact PNe by Miszalski et al. (2008). 

Modest variations in measured IAM stellar parameters from 
SHS data occur as a function of field position and the resulting 
variable PSF due to the effects of field rotation on the sky. This is 



7 http://www-wfau.roe.ac.uk/sss/halpha/ 



apparent when comparing point sources taken on the long 3 hour 
Ha exposures with the short 15 min SR exposures, and also from 
field vignetting, especially at large radii from the field centres. To 
better account for such effects requires that photometric compar- 
isons are performed over limited 1 -degree areas such that the true 
Ha emitters show a bona-fide enhanced Ha magnitude compared 
with the continuum SR magnitude. At the bright end of the mag- 
nitude distribution, severe saturation effects also come in to play, 
limiting stellar photometry to R of about ^12 in both the Ha and 
SR passbands. Further details are given in Paper I. 



7 CONCLUSIONS AND RECOMMENDATIONS 

We have established calibration factors (Table [TJ for 170 of the 233 
SHS fields (73 per cent), so these can be used directly to derive 
accurate integrated Ha fluxes for a wide range of objects includ- 
ing PNe, Hll regions, Wolf-Rayet shells and supernova remnants 
(refer to Figure [TT). Based on the significant number of calibrating 
objects we have used, we can demonstrate that our derived fluxes 
agree to within ±0. 1 dex with independent fluxes from the litera- 
ture over a remarkable range of 7 dex in integrated flux. There is 
a small systematic offset where SHS fluxes seem to be approxi- 
mately 10-15% fainter than the equivalent SHASSA fluxes. The 
origin of this offset is attributed to the differing filter transmission 
of the [N II] lines between SHASSA and the SHS, and the fact that 
the original analysis assumed that the diffuse emission from the 
ISM was solely Ha emission. Better agreement can be obtained by 
scaling up the SHS Ha fluxes by 0.05 dex. 
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Figure 12. The offset between SHS and literature Ha fluxes plotted as a 
function of Ha surface brightness (in erg cm -2 s -1 sr _1 ) derived from 
the literature flux. The twin effects of photographic saturation and flux 
over- subtraction occur for the brighter nebulae with a surface brightness, 
logS(Ha) ^ -4.0ergcm- 2 s- 1 . 



Determining fluxes for extended objects within the first- 
ranked SHS fields can be performed with confidence. This is to be 
expected since these fields exhibit a tight linear data-data relation 
with respective SHASSA images. For most of the second-ranked 
fields, the flux uncertainties are somewhat larger, though still us- 
able. There are currently 63 relatively featureless fields in Table [TJ 
for which we could not determine the CF. These fields either have 
high chemical fog levels or the fields have a poor dynamic range. 
This does not mean there is anything wrong per se with the expo- 
sures, just that the fields are mostly devoid of variable intensity Ha 
emission making it difficult to determine the CF. Such fields may 
still contain significant numbers of small localised emission neb- 
ulae such as PNe. In the meantime, we use the mean CF for all 
first-ranked fields of 12.0 counts px _1 R _1 as a default value for 
fields with indeterminate values. The Ha surface brightnesses and 
integrated fluxes so derived should be accurate to better than ±0.2 
dex. Many of these fields still could be used with more confidence 
if local CF values were established, even for fields affected by high 
chemical fog. 

The adopted value for the SF and CF and rank status of each of 
the 233 SHS fields are given in Table Q]and is the main legacy ma- 
terial from this work. The fields can be used not just for exploratory 
science and discovering new astronomical objects, but now for also 
deriving their integrated Ha fluxes and surface brightnesses. The 
SHS will remain useful for many years to come, and will be super- 
seded by the new VPHAS+ Ha survey (J. Drew et al., in prepara- 
tion) only in the region along the SGP within \b\ ^5°. Even there, 
the SHS will be the source of first-epoch properties for emission- 
line stars and other exotica for future studies on variability. 

In summary, the following notes and caveats should be con- 
sidered when using the SHS data to provide Ha intensities and in- 
tegrated flux measurements: 

1 . The pixel data has been treated by a filter flat-field as described 



in Paper I. This is a one-off correction that is assumed valid for 
all survey fields. 

2. We recommend that objects located more than 1.8° from the 
SHS field centre not be measured for flux if accuracy is a re- 
quirement. Excepting for the fields on the east, south and north- 
ern edges of the SHS zone, the generous overlap between field 
centres means that a neighbouring field should be available for 
flux measurement. Note that when full 0.67'Vpixel resolution ar- 
eas are obtained from the SHS webpage using the "Get Image" 
tool, the best data is automatically selected from the unvignetted 
central region of an available field. 

3. The effect of geocoronal Ha emission should be low. Paper I 
found from a random sample of six SHS fields that the Earth 
shadow height for the whole 3-hour observation and across the 
five-degree field of view was greater than 6000 km, so low-level 
geocoronal emission (<2 R) is not problematic, as most SGP 
fields covered by the SHS will contain significantly stronger dif- 
fuse emission. 

4. Telluric OH emission may contaminate the images at low flux 
levels, but this is also altitude dependent and can be safely ig- 
nored. 

5. The image data is non-linear at high intensities. Photographic 
saturation sets in at around 500 - 600 R, or approximately 6000 
- 7000 counts per pixel (for a typical CF value). Objects show- 
ing values above this, or showing 'negative' residuals after con- 
tinuum subtraction should be treated with caution, and the total 
fluxes can be considered to be firm lower limits. Consider point- 
source algorithms if the nebula is bright and compact. 

6. Brighter nebulae that are too extended for the application of 
point- source algorithms should be conspicuous enough to be 
detected in SHASSA. We recommend the use of this survey in 
such cases. 

7. Integrated fluxes accurate to ±0.1 - 0.15 dex, can be generally 
obtained across the survey fields, with the second-ranked fields 
have slightly larger systematic errors. 

8. Pending a more detailed analysis, a CF of 12.0 
± 3.0 counts px _1 R _1 can be used as a default value for 
the 63 fields with no formal determination (see Figure [7}. The 
integrated Ha fluxes for any objects derived from these fields 
should be accurate to better than ±0.2 dex. Alternatively, any 
known PNe or other small nebulae with an independent flux can 
be used to calibrate the whole field, but only if they are below 
the saturation limit. In such cases, the CF value for the field can 
be treated as a free parameter and varied until the measured 
integrated SHS flux matches the independent value. 

9. Emission from the red [Nil] lines often contributes signifi- 
cantly to the Ha+[N II] flux, especially for PNe and SNRs. This 
contribution needs to be carefully de-convolved from the total 
flux using the observed spectroscopic ratio for individual ob- 
jects, and following the procedure described in £14.31 

10. Owing to the different passbands of the SHASSA and SHS 
filters, the SHS Ha fluxes determined from Equation|2]should be 
increased by an additional factor; i.e. the logarithmic Ha fluxes 
should be brightened by 0.05 dex. 

11. The nature of photography and the scanning process produces 
a range of artefacts and spurious images on the SHS. The user 
should be aware of these, and can find more details in Paper I. 
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